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(57) Abstract: This invention is an audio system capable of identifying 
frequencies of feedback signals and filters these feedback signals. 
Frequency interpolation is utilized on a sampled frequency spectrum 
signal corresponding to a feedback signal, where the frequency 
interpolation allows the frequency of the feedback signal to be identified, 
especially where the frequency lies between samples of the frequency 
spectrum signal. The frequency identification allows a less intrusive 
filter, such as a notch filter, to be placed at the determined frequency 
of the feedback to eliminate the feedback signal. The placement of the 
notch filter reduces the effect on the audio signal provided by the audio 
system. The audio system may provide adaptive filtering of multiple 
feedback signals, using a single filter such as a notch filter. The single 
notch filter is configured to filter the more than one feedback signal, 
thereby freeing-up notch filters for filtering other feedback signals, or 
for reducing memory and processing requirements for the processor of 
the system. 
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AUDIO FEEDBACK PROCESSING SYSTEM 

5 

Background Of The Invention 

10 1. Related Applications. 

[01] This application claims priority to U.S. Provisional Patent Application 

Serial No. 60/363,994, entitled "Employing Narrow Bandwidth Notch Filters In 
Feedback Elimination," and filed on March 13, 2002, and is incorporated by 
reference. 

15 

2. Technical Field. 

[02] This invention generally relates to feedback in audio systems. More 

particularly, this invention relates to identifying a frequency of feedback and 
adaptive filtering of feedback signals in an audio system. 

20 

3. Related Art. 

[03] An audio system typically includes an input transducer (microphone), an 

amplifier, a microprocessor and an audio output (loudspeaker). The input transducer 
receives sound into the system, the amplifier amplifies the sound, the microprocessor 

25 performs signal processing, and the audio output (loudspeaker) provides sound to 
users of the system. Many audio systems allow for a duplex operation, where sound 
may be input to the microphone while audio is provided at the speaker. However, 
when the microphone receives a portion of the audio provided at the speaker as an 
input, an unstable, closed-loop system is created, resulting in audio feedback. 

30 [04] Audio feedback is manifested as one or more audio feedback signals at 

the speaker, where each feedback signal may be modeled as a sinusoidal signal (i.e. 
the feedback signal(s) exhibit characteristics of a sinusoidal signal). To eliminate a 

1 
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pabular feedback signal, the microprocessor converts (fie audio signal into a 
dtscrete (sampled) frequency spectrum representation, such as a Discrete Fourier 
Transform (DFT), Spectral Estimation, Filter Banks, or like representation The 
conversion of the audio signal to the sampled frequency specttum aHows for a 
genera, identification of the frequency of the feedback signal. The frequency samp!e 
havtng the greatest magnitude in the discrete frequency domain is selected as the 
frequency of the feedback signal. 

I0S1 A notch filter i s placed a. th. identified frequency of the feedback signal 
to ehminate t ba, pa rti cular feedback signal. However, because of computational and 
memory limitations of the microprocessor, the sampling resolution of the sampled 
frequency specfrum representation is limited. Thns, the selected frequency sample 
does no, provide an accurate estimate of an actual frequency orthe feedback signal 
Because the selected frequency sample is no, an accurate es,ima,e, a noteh filter is 
Utihzed ,ha, has a significantiy wider bandwidth and/or a greater cntidepth than what 
ts actually necessary for filtering the feedback signal. The wider bandwidth and/or 
greater cut-depth are necessary to ensure «ha, the feedback signal is ehminated from 
the otupu, sig„ aI . „ owever> me „ se rf a ^ ^ ^ ^ 

notch filter can degrade the audio quality of the sound at the speaker. 
|06| The computational and memor, fimitations of the micreprocessor limits 
*. t number of notch filtere mat may be used to eliminate audio feedback signa.s 
Where the number of feedback signals exceeds the number of notch filters available' 
some of the feedback signals camto. be eliminated by the system. The failure to 
ehmmate a, teas, some of the feedback sisals may require a system gain to be 
reduced, resulting in degraded system perfomtance. 
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Summary 

[07] This invention provides an audio system that identifies the frequency of a 

feedback signal using interpolative feedback identification. The interpolative 
5 feedback identification may be accomplished using frequency interpolation on a 
sampled frequency spectrum signal corresponding to a feedback signal. The 
feedback interpolation allows the frequency of the feedback signal to be identified, 
especially where the frequency of the feedback lies between samples of the 
frequency spectrum signal. The interpolation may include using samples of the 

10 sampled frequency spectrum signal to generate a unique quadratic (or higher order 
polynominal), which resembles the original main lobe of the feedback signal 
represented by the frequency spectrum signal. The polynominal may be constructed 
from the samples using polynominal interpolation, rational function interpolation, 
cubic spline interpolation, and the like. The peak of the polynominal and thus a 

15 representation/estimation of the actual frequency of the feedback signal may be 
determined, for example, by setting a derivative of the generated polynominal 
equation to zero. A narrowly tailored filter, such as a notch filter, may be placed at 
the determined frequency of the feedback to eliminate or reduce the feedback signal. 
The filter also reduces the effect on the audio signal quality provided by the audio 

20 system. 

[08] The audio system may adaptively filter multiple feedback signals using a 
single filter such as a notch filter. The adaptive filtering may include identifying 
frequencies of feedback in the audio signal, and determining which frequencies of 
feedback signals lie within a frequency window comprising adjoining samples of the 

25 sampled frequency spectrum. A filter, such as a notch filter is configured to filter out 
the frequencies identified as within a frequency range covered by the frequency 
window, thereby freeing-up notch filters for filtering other feedback signals, or for 
reducing memory and processing requirements for the microprocessor of the audio 
system. The frequency range covered by the frequency window may comprise any 

30 number of adjoining samples, and may be predetermined and/or configurable. 
Further, the frequency range covered by the frequency window may vary depending 
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on .he frequency band being examined, and/or ,be reso.ulion of ,be samp.ed 
frequency spectrum. 

I09| Olher systems, melhods, fealures and advances of Ihe invention „i„ be 
or w,„ beeome, apparen, ,o one wim in me art upon examina|ion ' 
fo< owing fi gures and delai|ed descrip ,. on J( , s imendKj tha( ^ such ^ 

systems, melhods, fealures and advantages be i„ cl uded within this description be 
wirtnn fte scope of the invemion, and be pro.ee.ed by the lowing elaima. ' 

Brief Description Of The Drawings 
110) The invention ean be heller onderslood wilt, referenoe lo Ihe Mowing 
drawings and deseriplion. The component in ,he figures are no, neeessan.y ,o sca,e 
emphasis inslead being p,aced upon iUuslraling Ihe principtes of ,he invemion" 
Moreover, in lhe flgures> , jke K[sr ^ • 

throughout the different views. 

HI J Figure, is . block diagram of m ^ ^ 

1 5 identification and reduction techniques. 

112] Figure 2 is a flow chart iterating operation of the audio system of 

F.gure 1 m identifying the frequency of a feedback signal. 

.13) Figure 3 is a graph illustrating a time-domain feedback signal 

20 I"! „ FiSUrC 4 " 3 *** i,IUStratin8 ^ DiSCrCte Time Fourier T ^orm of the 
feedback signal of Figure 3. 

.15] Figure 5 is a graph illustrating a time-domain window function 

116] Figure 6 is a graph illustrating a Discrete Time Fourier Transform of the 

time-domam window function of Figure 5. 

H7] Figure 7 is a graph illustrating the time-domain signal resulting from 



25 



, . - signal icsuiung ri 

multiplymg ,he feenbach signs! of Fignre 3 wilh Ihe window function of Pigure 5 

Figure 8 is a graph i U u sl ra«in g me Discrete Time Fourier Transform of me 
windowed feedback signal of Figure 7. 

U9] Figure 9 is a graph illustrating the Discrete Fourier Transform of the of 
the windowed feedback signal of Figure 7. 
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[20] Figure 10 illustrates an expansion of a portion of the graph of Figure 9, 
showing frequency bins which may be utilized in interpolating a frequency of a 
feedback signal. 

[21] Figure 11 is a graph comparing characteristics of prior art notch filters 

5 with a notch filter configured using interpolative feedback identification. 

[22] Figure 12 is another graph comparing characteristics of a prior art notch 

filter, with a notch filter configured using interpolative feedback identification. 
[23] Figure 13 is a flow chart illustrating operation of the audio system of 
Figure 1 for performing adaptive filtering. 
10 [24] Figure 14 is a graph illustrating a frequency window covering a specified 

frequency range for a time-domain signal, which may be utilized in performing 
adaptive filtering. 

[25] Figure 15 is a graph illustrating a frequency window covering a specified 

frequency range for a frequency-domain signal, which may be utilized in performing 
15 adaptive filtering. 

[26] Figure 16 is a graph illustrating characteristics for two notch filters for 

filtering corresponding feedback signals. 

[27] Figure 17 is a graph illustrating characteristics of a notch filter configured 

for adaptively filtering two feedback signals. 

20 Detailed Description of the Preferred Embodiments 

[28] Figure 1 is a block diagram of an audio system 100 having feedback 

identification and feedback reduction or elimination techniques. The audio system 
uses interpolative feedback identification and may adaptively filter multiple feedback 
signals using one notch filter. The interpolative feedback identification provides for 

25 a single estimate of the feedback frequency achieved from more than one sample of a 
discrete frequency spectrum representation of a feedback signal. The interpolative 
feedback identification may include utilizing frequency interpolation by generating a 
second degree or higher polynominal using one or more samples of the discrete 
frequency spectrum representation. An accurate representation of the actual 
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frequency of the feedback signal may be determined, for example, by setting a 
denvafve of the polynomial to zero. A filter, such as a notch filter, may be pl aced 
hi response to the interpolate feedback identification to reduce or eliminate the 
feedback signal with little or no effect on the audio signal quality provided by the 
5 aud,o system. The adaptive filtering involves configuring a filter, such as a notch 
filter, to eliminate multiple feedback signa.s, which allow other filters to reduce or 
ehmmate other feedback signals. The adaptive filtering may also, or in the 
alternative, reduce processor memory and/or computational requirements of the 
audio system. 

10 |29] The audio system 100 includes an audio input, i.e. a microphone 102 for 
receding an audio signal. The microphone 102 is coupled with a microprocessor 
104, wh,ch is capable of controlling operation of the audio system 100 The 
nucroprocessor 104 may perform any analog to digital conversions of audio signals 
receded and digital signal processing. The microprocessor 104 is further capable of 
15 performmg digital to analog conversions of audio provided by the audio system 100 
The microprocessor 104 is coupled with an amplifier 106 capable of amplifying an 
output audio signal. Amplifier 106 is coupled with a loudspeaker 108 for providing 
the output audio signal to a user of the audio system. While a particular 
configuration is shown, the audio system may have other configurations, including 
20 those with fewer or additional components. 

.30] Figure 2 is a flow chart of a method for identifying and reducing and/or 
removmg a feedback signal in an audio system. A time-domain audio signal s[n] 
from the microphone 104 is received 200 at microprocessor 104. Audio feedback 
may result when one or more portions of the audio provided from loudspeaker 108 is 
reccved at microphone 102, thereby causing an unstable, closed-loop system 
M.croprocessor 104 converts 202 the time-domain audio signal into a sampled 
frequency-domain signal \ S (k) . Microprocessor 104 may use windowing techniques 
such as Rectangular, Hamming, Bartlet, and like techniques to compute the 
frequency domain signal. The microprocessor 104 may then detect 204 the 
feedback. The detection of feedback may include performing frequency spectrum 



25 
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analysis such Discrete Fourier Transform (DFT), Spectral Estimation, Filter Banks, 
and like techniques. Samples of the frequency domains signal may be used in 
interpolating 206 to determine the frequency of the feedback signal, and the feedback 
signal may be filtered 208. Interpolating 206 and filtering 208 will be discussed 
5 further below with respect to Figure 10. 

[31] Figures 3-10 illustrate detecting of the feedback signal by microprocessor 

104. Figure 3 illustrates a time-domain feedback signal s[n]. Figure 4 illustrates a 
frequency domain signal resulting from converting the feedback signal s[n] 

to the frequency domain using, for example, the Discrete Time Fourier Transform 
10 (DTFT). Figure 5 illustrates a time-domain window function w[n]. Figure 6 
illustrates the DTFT (|^(^ y>r ]j ) of the window function w[n]. Figure 7 illustrates the 

product of the time-domain feedback signal s[n] with the time-domain window 
function w[n]. Figure 8 illustrates the windowed frequency domain signal |S(e /w )l 
centered about the frequency domain feedback signal |S , (e yV )|, resulting from taking 

15 the DTFT of the product of s[n] and w[n]. Figure 9 illustrates the sampled frequency 
domain signal \S[k]\ resulting from taking the DFT of the product of s[n] and w[n]. 
This is, for example, equivalent to sampling the windowed frequency domain 
feedback signal l^e 7 *)] of Figure 8 at equally spaced frequency intervals. Figure 10 
illustrates a portion of the sampled, windowed frequency domain signal \s[k\ of 

20 Figure 9, specifically showing a more detailed view around a main lobe of the 
feedback signal. The frequency spectrum signals illustrated in Figures 4, 6 and 8 are 
DTFT. The frequency spectrum signals illustrated in Figures 9 and 10 are DFTs. 
Other frequency spectrum analysis techniques may be utilized in converting the time- 
domain signal to the frequency domain, and analyzing the frequency domain signal. 

25 [32] In the flowchart of Figure 2, the interpolating 206 provides a single 

representation/estimation of a feedback frequency determined from multiple samples 
of the discrete frequency spectrum representation of the frequency signal. The 
interpolative feedback identification may be determined using frequency 
interpolation techniques, for example, as will be described with respect to the graph 
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of Figure 10, where each frequency sample defines a frequency bin. The notations 
used in Figure 1 0 are as follows: 

l 33 l Ktna* = The estimated frequency of the feedback signal. 

[34] B p = Peak (maximum) bin number. 

5 135] B p _, = Bin just below (in frequency) the peak bin number, 

f 36 J B P +, = Bin just above (in frequency) the peak bin number. 

I 37 J A es* ma ,e = Amplitude at the estimated frequency of the feedback. 

[38] A p = Amplitude of the peak bin. 

I 3 *] A p _, = Amplitude of the bin just below (in frequency) the peak bin. 

10 ! 4 °1 A p+t = Amplitude ofthebinjust above (in frequency) the peak bin. 

[41 J B^,, is the estimated frequency of the feedback signal which may be 
determined using the interpolation techniques described below. Ideally, the 
frequency B„, /ma , e will coincide with the actual frequency of the feedback signal. In 
any event, the frequency B estimale is typically a more accurate estimate of the actual 
frequency of the feedback signal than the frequency B p which is selected by systems 
of the prior art. 

[42] Interpolate feedback identification such as frequency interpolation 
provrdes a more accurate estimate of the actual frequency of feedback, and may be 
determined using samples of the DFT \ S [k) . Using the samples of the DFT signal 
a unique quadratic (or higher order polynomial) may be generated which 
resembles the original main lobe of the DTFT representing the feedback signal. A 
polynomial may be reconstructed from the sample points of the DFT \ S [k] . An 
interpolating polynomial for degree N-l is illustrated as a LaGrange polynomial by: 

25 Other mterpolating polynomial techniques may be used, including polynomial 
interpolate, rational function interpolation, cubic spline interpolation and the like. 
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[43] Applying the LaGrange polynomial equation to frequency interpolation 

(here, for a 2 nd order quadratic) yields a feedback frequency equationyfB) of: 

(B - B )(B - B j ) (B - B j )(B - B § ) (B - B , )(B - B ) 

f(B) = ~ — A p _j + — — A p + — A p+J 

(B^j -B p )(B p _j -B^j) (B p -B^HBp-B^,) (B ^ - B ^ )(B ^ - B p ) 

5 [44] A peak of the quadratic curve, and thus an estimate/representation of the 

frequency of the feedback signal may be determined by solving for a maximum of 
/(B). Solving for the maximum may be accomplished, for example, by taking the 
derivative of/(B), and setting the derivative to zero, yielding the estimated feedback 
frequency B es timate as: 
10 [45] 

, K,/ * ( B P + B p+> )<*p- B^)(B p -B p ^)(B p+l -B p _ I )(B p+l -B p )\ 



15 



B estimate 



| [A p *(B M -hB^KB^ -B p )(B M -B^XB^ -B M )(B p+1 -B p )j 

2 



[A p+1 * (B p _, + B p )(B^ - B p )(B p _ 1 - B p „ )(B p = B p _, )(B p - B p „ )j 

2 



[46] The pole of the quadratic curve provides a more accurate representation 

of the frequency of the feedback signal than the frequency B p of the peak bin alone. 
Where it is known that, prior to the interpolation, A p is greater than both A p+ i and A p . 
20 i, it may be determined that the interpolated polynomial has no minimum at this 
location, and only a maximum. Thus, taking the derivative of the interpolation 
polynomial and setting it to zero yields the maximum, and thus the more accurate 
representation of the frequency of the feedback signal than the frequency B p . 
However, where it is not known prior to the interpolation that A p is greater than both 
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lzz d m Ap : *r* necessary to detennine that the frequency at b — * * 

maximum and not a minimum of the quadratic equation. 

H7J To determine that the frequencv at Ft ;„ 

ic ircquency at B estjmatc , s a maximum (and not a 

m.„,mum) of the quadratic equation, a va,ue W may be compute. by the 
nncroprocessor ,04 using ,„e equation for /(*) abo ve, repenting ,„e amphnode 

compared w„h nte va,ues Ap*| aud A H , which are ampntudea of ,he feedback signa, 
^corresponding frequencies B, and ^. t0 ensure tha , A _ ^ ^ ^ 

'0 148] The inrerpoiating 206 of Figure 2 provides a more accurate estimate of 

Alter may be configured f„, fiftering 208 the feedback of ,he audio signa,. The finer 
may be a bandwidth finer. Other fine. may be used. Since a Cose estimate 
for me frequency of the feedback signa, has been identified using frequency 
■n eg ation me bandwidth noteh fi.ter may he configured (i.e., coefficient 
ca,cu,a.ed .herefore inCnding Quati.y Factor and/or g ain/cu,-dep,h, by .he 
nncroprocessor ,04 as a na^w bandwidth noteh fiher capab,e of fi ltering . oul lhe 
frequency of me feedback signa,. The microprocessor ,04 may a,so minimize a, 
■ea* one of a ba„d„id,h and a cu.-depU, of ,he noteh fi„er. The configured fiher 
may men be p,aced a, me frequency B B „ ra „ ,i. e . designed wi,h a center frequency of 
W). Such fihering may he emp.oyed utihzing fihering tec hni,„es such as Finite 
lmpu,se Response (FIR) and Infinite Impu.se Response (HR) technique,, or any 
o,her fi.tenng technique sufficient for fihering on, the feedback signa, as wou,d be 
apprectated by one skifieo in the art. Thus, identifying the frequency of the feedback 
stgna, usmg intensive feedback identification a„o„s for more accurate p.acemen, 
of the notch fi„er a, the frequencyof.be feedback signa,, and «hns is more accurate.y 
configured for fi,tering-ou. lhe feedback signal. 

149] Figure ,0 ifiushates an examp.e of intention by generating a 
po.ynom.a wh.ch mode,s .he origina, main ,obe of ,he frequency spechum, where 
■he ,n,e,o,a,.on is camed-o„, by sowing for a maximum of the po.ynomia, by 
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derivation. One skilled in the art would realize that any interpolation techniques may 
be utilized to identify the feedback frequency. For example, additional frequency 
bins may be interspaced between samples of the sample frequency domain signal 
shown in Figure 10, each interspaced bin having zero energy value. The sampled 
5 frequency domain signal may then be passed through a low pass filter resulting in an 
interpolated sampled spectrum. Using the interpolated sampled spectrum, one could 
identify a maximum of the filtered frequency spectrum to obtain a more accurate 
estimate of the feedback signal frequency. 

[50] Figures 1 1 and 12 illustrate graphs comparing characteristics of prior art 
10 notch filters with notch filters configured in accordance with interpolative feedback 
identification. The sampled frequency bin having a maximum amplitude B p in 
Figure 10, may correspond to 994 Hz in Figures 11 and 12. A more accurate 
representation of the frequency of the feedback signal, Estimate in Figure 1 0, may 
correspond to 1000 Hz in Figures 11 and 12. The sampled frequency bins and 
15 frequency of the feedback signal may have other frequencies. As shown at 
Figures 1 1 and 12, prior art feedback identification techniques result in a notch filter 
being configured to filter out frequencies at the maximum bin frequency 994 Hz, and 
thus must have an increased bandwidth as shown by line 1100 Figure 11, or 
increased cut-depth as shown by line 1200 of Figure 12, to ensure that the gain (G) of 
20 the filter at the actual frequency of the feedback is sufficient for filtering the 
feedback signal. 

1 51] In contrast, feedback identification techniques using interpolative 
feedback identification provide a more accurate representation (here about 1000 Hz) 
of the actual frequency of feedback. Accordingly, a notch filter having 

25 characteristics shown at 1105 and 1205 of Figures 11 and 12 may be placed at the 
more accurate estimate for the actual frequency of the feedback signal. Because the 
filter is more accurately placed, it may be more narrowly tailored (i.e. reduced 
bandwidth and/or cut-depth) while ensuring that the gain is sufficient at the 
frequency of the feedback signal to eliminate or reduce the feedback signal, and 

30 having little or no effect on the quality of the signal provided at the loudspeaker 108, 

11 

BNSDOCID: <WO 03079721A1_I_> 



WO 03/079721 



PCT/US03/07943 



10 



15 



20 



25 



30 



or in any event, less of an effect on the audio quality than notch filters configured 
using prior art feedback identification techniques. 

152] Figure 13 is a flow chart of a method for providing adaptive filtering of 
feedback in an audio system. Frequencies of a plurality of feedback signals are 
identified/estimated 1300 by the microprocessor 104. Such frequencies may be 
identified as described above using interpolate feedback identification, or in any 
other fashion. The microprocessor 104 determines 1302 whether the frequencies of 
feedback signals are within a frequency window covering a specified frequency 
range. The frequency range covered by the frequency window may be 
predetermined and/or configurable, and may vary depending on the frequency band 
being examined. The specified frequency range covered by the frequency window 
will be discussed further below with respect to Figures 14 and 15. 
153] The microprocessor 104 filters 1304 the feedback signal within the 
frequency range covered by the frequency window. The microprocessor 104 
configures a filter for filtering out any frequencies a feedback signal determines to be 
within the frequency range. The filter may be a notch filter or other type of filter. 
The microprocessor may determine filter coefficients such as quality factor, cut- 
depth and a center frequency for the filter. 

[54] Figure 14 is a graph illustrating a frequency window covering a specified 
frequency range for time-domain representations of feedback signals, which may be 
utilized in providing the adaptive filtering discussed above with respect to Figure 13. 
As shown in Figure 14, a frequency window represented generally at 1405 may 
cover a specified frequency range, for example, af. Where two feedback 
frequencies, for example feedback frequency fl and feedback frequency f2 lie within 
the frequency window 1405, it may be determined 1302 that adaptive filtering will 
be utilized to configure a single filter to filter out the feedback frequencies. 
[55] To determine whether the feedback frequencies lie within the frequency 
window 1405, a frequency differential Af may be determined between feedback 
frequencies, for example by subtracting one frequency from another. For example, 
as shown in Figure 14, Af may be determined by subtracting the frequency fl 
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representing a first frequency at which feedback is located from £2 representing a 
second frequency at which feedback is located. Where the value Af is less than af, 
and thus the frequency range covered by the frequency window 1405, it may be 
determined that the feedback located at frequencies fl and £2 may be adaptively filtered 
5 by a single filter. 

[56] A filter may be configured, for example by the microprocessor 104 at a 

center frequency fc within the frequency window 1405 having sufficient quality factor 
and/or cut-depth to filter out the feedback at the frequencies fl and £2. 
[57] Concurrently or subsequently, if a feedback signal is identified as being 

10 located at a frequency O, for example as shown in Figure 14, the microprocessor 104 
may determine whether the frequency differential Af between f3 and fc is less than the 
frequency range af covered by the frequency window 1405. Where it is determined 
that the newly calculated Af is less than af, the microprocessor 104 may determine that 
the feedback identified at f3 may be adaptively filtered utilizing the filter at fc, and thus 

15 reconfigure the filter centered at fc (i.e., reconfigure the quality factor, cut-depth and/or 
fc) to filter out the feedback identified at the frequencies fl, £2 and £3. 
[58] Alternatively, instead of determining the frequency differential between f3 

and fc, the microprocessor 104 may instead determine a frequency differential Af 
between O and fl for comparing with the frequency range af of the frequency window 

20 1405 in determining whether the feedback frequencies fl , f2 and f3 may be adaptively 
filtered by a single filter. As additional feedback frequencies are concurrently and/or 
subsequently identified, the microprocessor 104 may determine whether to employ 
additional filters, or to utilize existing filters to cover the concurrently or subsequently 
identified frequencies of feedback. 

25 [59] In addition, the microprocessor 104 may further utilize algorithms that may 
minimize the number of filters necessary to filter out the identified feedback 
frequencies. In Figure 14, the frequency of the feedback frequency fl may be 1000Hz, 
where the feedback frequency £2 may be 1012 Hz and the feedback frequency £3 may 
be 1024 Hz. The specified frequency range af of the frequency window 1405 may be 

30 any value, for example, 6 Hz, 12 Hz, 20 Hz, 100 Hz or any other value. The specified 

13 

BNSDOCID: <WO 03079721A1 J_> 



WO 03/079721 



PCT/US03/07943 



10 



frequency range af may vary across the frequency spectrum, as a function of the 
frequency of the particular feedback frequencies being examined. For example, the 
frequency range af may increase logarithmically as the particular frequency being 
examined for feedback increases. Thus, at lower frequencies, af may have a smaller 
value than af at higher frequencies. In addition, the value of af defining the frequency 
window 1405 may be configurable by a user of the system 100. 

160] The graph of Figure 14 describes how the determining 1302 may be 
accomplished for feedback signals represented in the time-domain. The determining 
1310 may similarly be carried-out for identified feedback signals in the frequency 
domain, for example as described with respect to the graph of Figure 15. 
I61J Figure 15 is a graph illustrating a frequency window covering a specified 
frequency range for frequency domain representations of feedback signals, which may 
be utilized for the adaptive filtration discussed above. A frequency window 1505 is 
shown, covering a specified frequency range represented by a particular number of 
frequency bins (i.e., frequency samples) aB. To determine 1302 whether the 
feedback frequencies lie within the frequency window 1505, a frequency differential 
represented here as a number of frequency bins, AB, may be determined between 
feedback frequency bins, for example by subtracting one feedback frequency bin 
from another. As shown in Figure 15, AB may be determined by subtracting the 
frequency bin# B328 representing a first frequency at which feedback is located from 
the frequency bin# B326 representing a second frequency at which feedback is 
located. Where the value AB is less than aB, and thus the frequency range covered by 
the frequency window ,505, it may be determined that the feedback located at 
frequency bins B328 and B326 may be adaptively filtered by a single filter. 
[62] A filter may be configured, for example by the microprocessor 104 at a 
center frequency fc within the frequency window 1 505 having sufficient quality factor 
and/or cut-depth to filter out the feedback at the frequency bins B326 and B328. 
[63] Concurrently or subsequently, if a feedback signal is identified as being 
located at a frequency bin #B333, for example as shown in Figure 15, the 
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microprocessor 104 may determine whether the frequency differential AB between the 
frequency bin #B333 and fc is less than the specified frequency range otB covered by 
the frequency window 1505. Where it is determined that the newly calculated AB is 
less than aB, the microprocessor 104 may determine that the feedback identified at 
5 frequency bin #B333 may be adaptively filtered utilizing the filter at fc. The 
microprocessor 104 may reconfigure the filter centered at a center frequency fc (i.e., 
reconfigure the quality factor, cut-depth and/or fc) to filter out the feedback identified at 
the frequencies represented by frequency bins 326, 328 and 333. In Figure 15, the 
center frequency fc is shown, by example, at bin #B327. 

10 [64] Similar to as discussed above with respect to Figure 14, instead of 

determining the frequency differential between bin #B333 and fc, the microprocessor 
104 may instead determine a frequency differential AB between bins B333 and B326. 
This frequency differential AB may be compared with the frequency range aB of the 
frequency window 1505 to determine whether the feedback frequencies represented at 

1 5 bins B326, B328 and B333 may be adaptively filtered by a single filter. As additional 
feedback frequencies are concurrently and/or subsequently identified, the 
microprocessor 104 may determine whether to employ additional filters, or to utilize 
existing filters to cover the concurrently or subsequently identified frequencies of 
feedback. 

20 [65] Additionally, and as discussed above, the microprocessor 104 may further 

utilize algorithms that may minimize the number of filters necessary to filter out the 
identified feedback frequencies. The specified frequency range aB of the frequency 
window 1 505 is shown in Figure 15 as being 3 frequency bins, where the bin #326 may 
represent a frequency sample at 1000Hz, and spacing between frequency samples/bins 

25 may be approximately 6 Hz. However, similar to as discussed above with respect to 
Figure 14, it will be appreciated by one skilled in the art that otB may be any number of 
frequency bins, for example 2, 3, 5 or 10 frequency bins, and that the frequency 
differential represented by aB may vary as a function of the feedback frequencies being 
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examined. In addition, the value of aB rfpfinino *u*> • ^ 

vaiue oi an defining the frequency window 1505 may be 

configurable by a user of the system 1 00. 

[661 Figure 16 illustrates a graph showing characteristics of adjacently placed 
notch filters that may benefit from the adaptive filtering discussed herein. Feedback 
has been identified at frequencies of fl equal to about 1000 and Q equal to about 
1 012 Hz. To eliminate the feedback identified at these frequencies, notch filters may 
be utihzed having the characteristics 1600 and 1605. The characteristics 1600 
■nclude a Quality Factor equal to about 128 and a cut-depth equal to about -6dB to 
eliminate or reduce the feedback. The characteristics 1605 include a Qualify Factor 
equal to about 128 and a cut-depth equal to about -6dB to eliminate or reduce the 
feedback. However, in utilizing adaptive filtering, microprocessor 104 is capable of 
determining that the frequency differential Af between feedback frequencies at 
frequencies fl and f2 are within a frequency range af defining a frequency window 
where af may be 15 Hz. Microprocessor 104 may configure a single notch filter to 
filter out the feedback from both identified feedback frequencies. 
167] In Figure 17, characteristics of a notch filter configured by the 
microprocessor 104 is shown at 1700. The characteristics indicate a notch filter 
designed for a center frequency fc of about 1006 Hz and having a Quality Factor of 
equal to about 45, and a cut-depth equal to about -6dB. The notch filter is placed 
between the two identified frequencies, here fl at about 1000 Hz and f2 at about 
1012 Hz, to filter out the feedback signal frequencies. The notch filter may be placed 
0.e. designed with a center frequency) at a midpoint of the frequencies of identified 
feedback, here about 1 006 Hz. The notch filter may be placed at any other frequency 
between the identified feedback frequencies, or within the frequency window being 
examined (not shown), sufficient for filtering out the identified feedback Where 
more than two frequencies of feedback signals are determined to fall within the 
frequency range af, an average frequency may be calculated for the determined 
frequencies of feedback, where the filter is placed at the average frequency 
Alternatively, a midpoint frequency between the greatest and lowest frequencies 
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determined to be within the frequency range <xf defining the frequency window may 
be selected for placement of the notch filter. 

[68] Thus, instead of requiring two or more notch filters to filter out multiple 

feedback signals within the frequency window defined by the frequency range af, a 

5 single notch filter may be utilized. Hence, the other notch filter(s) available in the 
audio system may be used to eliminate or reduce feedback at other frequencies. 
Rather than having additional notch filters, reducing the number of notch filters for 
filtering feedback signals may reduce the memory and/or processing requirements of 
microprocessor 104. The filtering may be accomplished as software executed on the 

10 microprocessor 104. 

[69] Further, multiple sets of frequencies of feedback signals may be identified 

by the microprocessor 104, where the microprocessor 104 configures a notch filter to 
filter the feedback signals corresponding to each set of feedback frequencies. 
[70] The audio system 100 discussed above may be utilized in cellular 

15 telephones, public address systems, speakerphones having duplex operation, or any 
other audio system that may suffer from feedback. The microphone 102 may be any 
input transducer sufficient for receiving audio into the audio system 100. The 
microprocessor 104 may be any microprocessor capable of performing the 
functionality/processing, including converting time-domain signals to sampled 

20 frequency domain signals. Further, although not shown, the microprocessor 1 04 may 
include, or may be coupled with, an external storage media such as computer 
memory that may include computer programming, executable on the microprocessor 
104, for carrying out one or more of the functionalities described herein. The storage 
. medium may be magnetic, optical or any other storage media capable of providing 

25 programming for the microprocessor 1 04. 

[71] The loudspeaker 108 may be any speaker capable of providing the output 

audio from the audio system 100. Alternatively, hardware components not shown 
may be coupled with the microprocessor 104 for performing the sampled frequency 
domain conversion where the microprocessor 104 does not possess such 

30 functionality. The filtering may be accomplished using software, hardware or a 
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cognation, and need not be limited to notch filtering techniques. The software 
may be executable on a microprocessor such as performing digital signal processing 
or the like. The hardware may be coupled with the microprocessor 104, which may 
configure the hardware to achieve desired processing and/or filtering characteristics 
(721 In addition, the values illustrated and discussed in relation to the Figures 
are exemplary, and are not limitations on the feedback identification and elimination 
or reduction system. Further, the value for the frequency range af with respect to 
adapt.ve filtering may be any value while achieving at least some of the advantages 
d.scussed herein. The frequency range affaB may be increased (made larger) to 
reduce the number of filters required to eliminate feedback. A lower number of 
filters may be desired where the number of feedback signals outnumber the number 
of filters available for filtering feedback, or where a processor performing the 
filtenng has limited memory and/or processing capabilities. The frequency window 
defined by the frequency range af/oB may be detennined based on considerations 
w.thm the particular audio system utilized, and may be user-configurable Such 
^derations may include selection of a frequency range which allows frequencies 
of feedback signals to be combined without unduly affecting the audio quality 
prov.ded by the audio system. However, different audio systems have varying 
reqmrements as to the audio quality provided thereby. For example, a public address 
system may have less stringent audio quality requirements than an audio system that 
may be used in a concert hall or the like. A larger frequency range value <xf7aB may 
be desired for the former than for the latter to account for desired audio quality. 
(73] Further, one skilled would realize that various techniques may be 
employed in identifying which frequencies of feedback within the frequency range 
aCccB. Further, the microprocessor may utilize various techniques in grouping 
.dent.fied feedback signal sets which are each to be filtered by a single filter, where 
the techmque may minimize the number of filters required for filtering the identified 
feedback signals. 

|74| The audio system 100 may perfotm bom intetpolative feedback 
tdentincation in identifying freouencies of feedback ^ ^ 
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for configuring a filter-to-filter out multiple frequencies of feedback signals. The 
audio system 100 need not perform the feedback identification using interpolative 
feedback identification and/or the adaptive filtering. Rather, the audio system 100 
may be utilized in identifying the frequencies of feedback using interpolative 
5 feedback identification while being coupled with additional hardware or 
microprocessing capabilities which are utilized in eliminating or reducing the 
identified frequencies of feedback. The hardware may include adaptive filtering. 
Further, the audio system 100 may perform adaptive filtering using frequencies of 
feedback identified by external hardware or a processing functionality (which may or 
10 may not include feedback frequencies identified using the interpolative feedback 
identification). 

[75] The illustrations have been discussed with reference to functional blocks 

identified as modules and components that are not intended to represent discrete 
structures and may be combined or further sub-divided. In addition, while various 
1 5 embodiments of the invention have been described, it will be apparent to those of 
ordinary skill in the art that other embodiments and implementations are possible that 
are within the scope of this invention. Accordingly, the invention is not restricted 
except in light of the attached claims and their equivalents. 



19 

BNSDOCID: <WO 03079721 A 1 J _> 



WO 03/079721 



PCT/US03/07943 



25 



Claims 

What is claimed is: 

1 • A method for processing an audio feedback signal, comprising: 
receiving an audio signal; 

5 applying an interpolate feedback identification to a feedback signal within 

the audio signal; and 

identifying a frequency of the feedback signal in response to the interpolate 
feedback identification. 

10 2. The method of Claim 1 , where the audio signal is a time-domain audio signal, 
and further comprising: 

windowing the time-domain signal; 

converting the time-domain signal to a frequency domain signal; and 
applying the interpolate feedback identification to identify the frequency of 
the feedback signal between samples of the sampled frequency domain signal. 

3. The method of Claim 2, further comprising generating a polynomial 
responsive to the samples of the frequency domain signal, and applying interpolate 
feedback identification responsive to the polynomial. 

4. The method of Claim 3, further comprising solving for a maximum of the 
polynomial to identify the feedback frequency. 

5- The method of Claim 4, further comprising verifying that the identified 
frequency has a greater magnitude than a magnitude of the frequency samples of the 
windowed frequency domain signal. 
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6. The method of Claim 2, further comprising: 

calculating coefficients for a notch filter responsive to the identified frequency 
30 of the feedback signal; and 

filtering out the identified frequency of the feedback signal using the notch 

filter. 
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7. The method of Claim 6, further comprising minimizing at least one of a 
bandwidth and a cut-depth of the notch filter responsive to the identified frequency of 
the feedback signal. 

5 

8. The method of Claim 1, further comprising filtering the audio signal to 
eliminate the identified feedback responsive to the interpolative feedback 
identification. 

1 0 9. The method system of Claim 1, where there are a plurality of feedback signals 
within the audio signal, and further comprising identifying a plurality of frequencies 
corresponding to the plurality of feedback signals responsive to the interpolative 
feedback identification. 

15 10. The method of Claim 9, further comprising filtering the audio signal to 
eliminate the plurality of frequencies of identified feedback responsive to the 
interpolative feedback identification. 

11. The method of Claim 9, further comprising adaptively filtering at least two 
20 adjacent identified frequencies of feedback signals responsive to the interpolative 

feedback identification with a single filter by configuring the filter to filter the at least 
two adjacent identified frequencies. 

12. The method of Claim 11, where the filter is a notch filter, and further 
25 comprising determining coefficients for the notch filter to minimize at least one of a 

bandwidth and a cut-depth of the notch filter responsive to the identified frequencies 
of the feedback signal. 

13. The method of Claim 1 1 , further comprising selecting the at least two adjacent 
30 identified frequencies of feedback within a specified frequency range. 
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14. The method of Claim 13, where the specified frequency range is variable 
responsive to a feedback frequency being examined. 

1 5. The method of Claim 1 1, further comprising: 

determining a minimum identified frequency and a maximum identified 
frequency which will be filtered; 

selecting an interior frequency between the minimum identified frequency and 
the maximum identified frequency; and 

configuring the filter at the selected interior frequency. 

1 6. An audio system comprising: 

an audio signal port for receiving an audio signal; 

a processor coupled with the signal port for applying interpolate feedback 
identification to a feedback signal within the audio signal, and identifying a frequency 
of a feedback signal within the audio signal responsive to the interpolative feedback 
identification. 
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17. The audio system of Claim 16, where the audio signal is a time-domain audio 
signal, and further comprising: 

the processor being capable of windowing the feedback signal, converting the 
windowed feedback signal to a frequency domain signal, and identifying the 
frequency of the feedback signal between samples of the sampled frequency domain 
signal responsive to the interpolative feedback identification. 



25 18. 



The audio system of Claim 17, where the processor identifies the frequency of 
the feedback signal between samples using the interpolative feedback identification 
by generating a polynomial responsive to the samples of the frequency domain signal. 

1 9. The audio system of Claim 1 8, where the processor identifies the frequency of 
the feedback signal responsive to the interpolative feedback identification by solving 
for a maximum of the generated polynomial. 
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20. The audio system of Claim 19, where the processor identifies the frequency of 
the feedback signal responsive to the interpolative feedback identification by 
verifying that the identified frequency has a greater energy magnitude than an energy 
magnitude of the frequency samples of the frequency domain signal. 

5 

21. The audio system of Claim 17, further comprising the processor calculating 
coefficients for a notch filter responsive to the identified frequency of the feedback 
signal. 

10 22. The audio system of Claim 21, where the processor calculates coefficients for 
the notch filter by minimizing at least one of a bandwidth and a cut-depth of the notch 
filter responsive to the identified frequency of the feedback signal. 

23. The audio system of Claim 16, further comprising a filter coupled with the 
15 processor and the audio signal port, for filtering the audio signal to reduce the 

feedback at the identified feedback frequency responsive to the interpolative feedback 
identification. 

24. The audio system of Claim 23, where the filter is a notch filter tuned to the 
20 identified frequency of the feedback signal. 

25. The audio system of Claim 16, where there are a plurality of feedback signals 
within the audio signal, and the processor applies interpolative feedback identification 
to the plurality of feedback signals, and identifies a plurality of frequencies 

25 corresponding to the plurality of feedback signals responsive to the interpolative 
feedback identification. 

26. The audio system of Claim 25, further comprising a plurality of filters coupled 
with the processor and the audio signal port, for filtering the audio signal to eliminate 

30 the identified feedback responsive to the interpolative audio. 
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27. The audio system of Claim 25, further comprising a filter coupled with the 
processor and the audio signal port, where the processor adaptively filters at least two 
adjacent identified frequencies of feedback signals by configuring a filter to filter the 
at least two adjacent identified frequencies, responsive to the inteipolative audio. 

28. The audio system of Claim 27, where the filter is a notch filter, and the 
processor configures the notch filter by determining coefficients for the notch filter to 
minimize at least one of a bandwidth and a cut-depth of the notch filter responsive to 
the identified frequencies of the feedback signal. 

29. The audio system of Claim 27, where the processor adaptively filters the at 
least two adjacent identified frequencies by selecting the at least two adjacent 
identified frequencies of feedback within a specified frequency range. 

30. The audio system of Claim 29, where the specified frequency range is variable 
responsive to a feedback frequency being examined. 

31 . The audio system of Claim 27, where the processor adaptively filters the at 
least two adjacent identified frequencies of feedback signals by determining a 
minimum identified frequency and a maximum identified frequency which will be 
filtered by the filter, selecting an interior frequency between the minimum identified 
frequency and the maximum identified frequency, and configuring the filter for 
filtering at the selected interior frequency. 

25 32. The audio system of Claim 27, where the processor includes a storage media 
programmed for applying the interpolative feedback identification to the feedback 
signal and identifying the frequency of the feedback signal responsive to the 
interpolative feedback identification. 
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33. A method of processing audio feedback, comprising: 

receiving an audio signal including multiple feedback signals; 

identifying a plurality of feedback frequencies, each feedback frequency 

corresponding to one of the feedback signals; 
5 determining whether at least two feedback frequencies of the plurality of 

feedback frequencies lie within a specified frequency range; and 

configuring a filter responsive to the determination to filter out the at least two 

determined feedback frequencies. 

10 34. The method of Claim 33, further comprising configuring a notch filter to filter 
out the at least two determined frequencies. 

35. The method of Claim 33, further comprising determining at least one of a 
bandwidth, cut-depth and center frequency of the filter. 

15 

36. The method of Claim 35, further comprising selecting a center frequency for 
the filter which is an average of the at least two feedback frequencies. 

37. The method of Claim 33 further comprising determining whether at least two 
20 adjacent feedback frequencies lie within the specified frequency range. 

38. The method of Claim 33, further comprising determining whether only two 
feedback frequencies lie within the specified frequency range. 

25 39. The method of Claim 33 where the specified frequency range is variable 
responsive to a feedback frequency being examined. 

40. The method of Claim 33, further comprising applying interpolative feedback 
identification to at least one of the feedback signals, and identifying at least one of the 
30 feedback frequencies responsive to the interpolative feedback identification. 
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41 - A storage media for use on a processor of an audio system, comprising: 

a memory portion programmed for allowing receipt of an audio signal 
applying an interpolate feedback identification to a feedback signal within the audio 
signal, and identifying a frequency of the feedback signal responsive to the 
interpolative feedback identification. 

42. The storage media of Claim 41, where the audio signal is a time-domain audio 
signal, and further comprising the memory portion being programmed for: 
windowing the time-domain signal; 

converting the time-domain signal to a frequency domain signal; and 
identifying the frequency of the feedback signal between samples of the 
sampled frequency domain signal responsive to the interpolative feedback 
identification. 

15 43. The storage media of Claim 42, further comprising the memory portion being 
programmed for generating a polynomial responsive to the sables of the frequency 
domain signal. 
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44 



The storage media of Claim 43, further comprising the memory portion being 
programmed for identifying the frequency of the feedback signal by solving for a 
maximum of the polynomial. 



45. 



The storage media of Claim 41, further comprising the memory portion being 
programmed for filtering the identified frequency of the feedback signal by 
calculating coefficients for the filter responsive to the identified frequency of the 
feedback signal. 



46. 



The storage media of Claim 45 where the filtering approximates a notch filter 
and further comprising the memory portion being programmed for minimizing at least' 
30 one of a bandwidth and a cut-depth of the approximated notch filter responsive to the 
identified frequency of the feedback signal. 
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47. A storage media for use on a processor of an audio system, comprising: 

a memory portion programmed for receiving an audio signal including 
multiple feedback frequencies, 

identifying a plurality of feedback frequencies, each feedback frequency 
5 corresponding to one of the feedback signal, 

determining whether at least two feedback frequencies of the plurality of 
feedback frequencies lie within a specified frequency range, and 

configuring a filter to filter out the at least two determined feedback 
frequencies responsive to the determination. 

10 

48. The storage media of Claim 47, further comprising the memory portion being 
programmed for configuring a notch filter to filter out the at least two determined 
frequencies. 

15 49. The storage media of Claim 48, where the notch filter is an approximated 
notch filter, and fiirther comprising the memory portion being programmed for 
approximating the notch filter. 

50. The storage media of Claim 47, where the memory portion is programmed for 
20 applying interpolative feedback identification on at least one feedback signal, and 
identifying at least one of the feedback frequencies responsive to the interpolative 
feedback identification. 
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